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Abstract. 

Starting from the Argonne V18 nucleon-nucleon interaction and using the Unitary Correlation 
Operator Method, a correlated interaction vucom has been constructed, which is suitable for cal¬ 
culations within restricted Hilbert spaces. In this work we employ the Vucom in Hartree-Fock, 
perturbation-theory and RPA calculations and we study the ground-state properties of various 
closed-shell nuclei, as well as some excited states. The present calculations provide also important 
feedback for the optimization of the Vucom and valuable information on its properties. The above 
scheme offers the prospect of ah initio calculations in nuclei, regardless of their mass number. It 
can be used in conjunction with other realistic NN interactions as well, and with various many-body 
methods (Second RPA, QRPA, Shell Model, etc.). 


The Unitary Correlation Operator Method (UCOM) provides a novel seheme for ear- 
rying out nuelear strueture ealeulations starting from realistie nueleon-nueleon (NN) in- 
teraetions iHllii . The major short-range eorrelations, indueed by the strong repulsive 
eore and the tensor part of the NN potential, are deseribed by a state-independent unitary 
eorrelation operator C. This ean be used to introduee eorrelations into an uneorrelated 
A—body state or, alternatively, to perform a similarity tranformation of an operator of 
interest. Applied to a realistie NN interaetion, in partieular, the method produees a "eor- 
related" interaetion, vucom, whieh ean be used as a universal effeetive interaetion, for 
ealeulations within simple Hilbert spaees. The same transformation ean then be applied 
to any other operator under study, as is needed for a eonsistent UCOM treatment. 

The utilization of the UCOM involves a eluster expansion of the eorrelated operators 
and, eurrently, a truneation at the two-body (2B) level. The latter is justified by the 
short range of the eorrelations treated by the method. The aim is to treat explieitly only 
the state-independent short-range eorrelations (SRC); long-range eorrelations (LRC) 
should be deseribed by the model spaee. The eorrelation funetions are determined 
for eaeh {S,T) ehannel by minimizing the energy of the two-nueleon system. Three- 
body interaetions are eurrently not ineluded in the formulation; one way to aeeount for 
these is by adding a simple phenomenologieal non-loeal 2B eorreetion to the eorrelated 
Hamiltonian. The introduetion of sueh a eorreetion (the same for all nuelei) has allowed 
the UCOM to sueeessfully deseribe properties of nuelei up to mass numbers A 60, in 
the framework of variational ealeulations within Fermionie Moleeular Dynamies |13]. 

In this work we study nuelear strueture and response, based on realistie interaetions, 
without restrieting ourselves to light-to-medium systems. This is made possible by 
employing the VucoM in Hartree-Foek (HF)- and RPA-based models. The 2B matrix 



FIGURE 1. Binding energy per nucleon, for the indicated nuclides, in HF (blue dots) and HF+PT (red 
squares) and experimental (bars). 


elements of the vycoM (without correction terms), Coulomb interaction and (intrinsic) 
kinetic energy are calculated in a harmonic-oscillator (HO) basis. These are used as 
input to subsequent HF and RPA calculations, in configuration space. Only spherical, 
closed-shell nuclei have been considered so far. The following methods have been used: 

HF - a spherical-HF method, to estimate at "zeroth order" the nuclear ground state (gs) 
properties. The HF single-particle states are expanded in the HO basis. 

HF-i-PT - second order perturbation theory (PT) is performed on the HF basis to obtain 
a correction to the gs energy. 

HF-i-RPA - a self-consistent model, which allows us to estimate a correction to the gs 
energy due to LRC, and in addition to study collective excitations. 

ERPA - an extended RPA |@], which is built on top of the true RPA gs and involves 
an iterative solution of the RPA equations. Corrected single-particle energies and 
occupation numbers can be obtained, as well as excitation properties. 

The PT and (E)RPA allow us to account for ERC, in addition to the SRC introduced 
via the UCOM. The optimal separation of the two types of correlations remains an 
important task. It is expressed primarily by a constraint on the range of the tensor 
correlation functions, imposed during their parameterization. By varying this range 
(more precisely, the "correlation volume" !■§ l^), a family of correlators and respective 
correlated interactions are obtained. Our present results provide important feedback for 
the optimization of the vucoM and valuable information on its properties. 

The maximum HO-energy and angular-momentum quantum numbers used here are 
Amax = 2n-\-l= 12 and ^max = 10, respectively, providing a satisfact(^ degree of 
convergence. We use the vucom parameterizations discussed in Ref. [ 0 ], based on 
the Argonne VI8 interaction. The results presented in Eigs. 1 and 2 were obtained 

with = 0.09 fm^. This value best reproduces, in exact no-core shell model 

calculations, the binding energy of the light systems ^He, ^H fl. 

Ground state properties. Binding is achieved already at the HE level. The inclusion 
of ERC corrections to the gs energy via PT brings the gs energy very close to the 
experimental data, as shown in Pig. 1. The RPA correlations result in overbinding |@], 
mainly due to the double-counting of second-order corrections inherent in the model. 
In general, larger-/^ tensor correlators provide stronger binding both at the HP and the 
RPA level. The HP single-particle levels (not shown) are too sparse |@]. The omission of 






FIGURE 2. Isoscalar monopole resonance of ^**Ca, in HF+RPA (blue lines) and ERPA (red lines). An 
arrow indicates the experimental centroid. 


a three-body interaction and of LRC are responsible for this effect. Occupation numbers 
Hi = (aj Qi) (in standard notation) have been calculated within ERPA for and "^®Ca. In 
principle, correlated operators should be used. The small depletion of the Fermi 

sea that we obtain using uncorrelated operators reflects the effect of the gs LRC 0] ■ 

Collective excitations. Our RPA results on the isoscalar (IS) giant monopole reso¬ 
nance (ISGMR), for various medium and heavy nuclei, are in good agreement with the 
experimental data. An example is shown in Fig. 2, where we see also that the gs cor¬ 
relations taken into account in ERPA have a relatively small effect. In general, lower 
7^ values result in lower ISGMR energies 0. The isovector (IV) dipole strength (not 
shown) is distributed at energies which are too high compared with experiment; ERPA 
is not able to correct for this result. Inclusion of 2p2h configurations within Second RPA 
is expected to bring the IVD strength to lower energies. 

In conclusion, the performance of the vucoM is very encouraging. Certain aspects 
of the model (e.g., optimal tensor correlators, three-body effects) are still open to im¬ 
provement. The above scheme can be used in conjunction with other realistic (local or 
non-local) NN interactions, as well as with various many-body methods (Second RPA, 
QRPA, Shell Model, etc), and offers the prospect of ab initio calculations across the 
nuclear chart. 
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